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We examine the static and dynamic stability of the solutions of the Gross-Pitaevskii equation and 
demonstrate the intimate connection between them. All salient features related to dynamic stability 
are reflected systematically in static properties. We find, for example, the obvious result that static 
stability always implies dynamic stability and present a simple explanation of the fact that dynamic 
stability can exist even in the presence of static instability. 
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INTRODUCTION 

The question of excitations is central to the study of 
cold atoms. Numerous experimental and theoretical in- 
vestigations have been devoted to the study of a variety 
of collective and elementary excitations in these gases in- 
cluding vortex states, solitary waves, and normal modes. 
When determining the properties of an excited state, it 
is natural to consider two kinds of possible instabilities 
- static and dynamic. In the former case, one wishes 
to determine whether a state which extremizes the en- 
ergy is a genuine local minimum of the energy subject to 
certain physically motivated constraints. Since the wave 
function extremizes the energy, infinitesimal perturba- 
tions will make no first-order change of the energy. The 
state will be stable provided that an arbitrary infinites- 
imal variation of the wave function necessarily increases 
the energy to second order. 

For the problem of dynamic stability, one conven- 
tionally considers the temporal evolution of arbitrary 
infinitesimal perturbations to the wave function using 
approximate linearized equations (e.g., the Bogoliubov 
equations). While the full time-dependent equations con- 
serve probability, the linearized equations do not, and 
the resulting eigenvalue problem is non-Hermitean. If the 
corresponding eigenvalues are real, the system is dynami- 
cally stable. In this case small-amplitude time-dependent 
perturbations of the solution lead to bounded motion 
about the original equilibrium position. Instability is 
indicated by complex eigenvalues, and the associated ex- 
ponential growth of small-amplitude perturbations drives 
the system to a new state beyond the scope of the lin- 
earized equations. 

Here, we will explore the intimate connection between 
these two seemingly-different criteria for stability. Two 
main points will emerge. First, static stability implies 
dynamic stability, but dynamic stability can exist even 
in the absence of static stability. Second, the transition 
from dynamic stability to instability is reflected in the 
features of the corresponding problem of static stability 
(subject to appropriate constraints). Here, inspired by 



the experiment of Ref. 'I'l and by the numerical studies 
of Refs. we will focus on the specific problem of 

the stability (static and dynamic) of a doubly-quantized 
vortex state. Our study is, however, more general, and 
our results can be applied to the study of the stability of 
any problem described by the non-linear Gross-Pitaevskii 
equation. One remarkable observation obtained from nu- 
merical simulations of this problem 2, 3, 4] is that the 
system alternates between regions of dynamical stability 
and instability as the strength of the interatomic inter- 
action increases. Since some (and probably all) regions 
of dynamical stability coincide with regions of static in- 
stability, it is useful to seek a simple description of this 
surprising phenomenon. 

In the following we first examine the questions of static 
and dynamic stability separately. We will then demon- 
strate the connection between them. 



MODEL 

We consider N atoms subject to a spherically- 
symmetric single-particle Hamiltonian, hi, interacting 
through a short-range effective interaction, 

N N 

H^Y.^., + Uo J2 '^(i-^ - r,)/2. (1) 

Here Uq = inh'^a/M is the strength of the effective two- 
body interaction, a is the scattering length for elastic 
atom-atom collisions, and M is the atomic mass. 

For simplicity we assume strong confinement along the 
z axis, which is also chosen as the axis of rotation. This 
assumption implies that the cloud is in its lowest state of 
motion in the z direction, and the problem thus becomes 
two-dimensional with higher degrees of freedom along the 
z axis frozen out. We also assume that hi is rotationally 
symmetric about the z axis (e.g., a harmonic oscillator 
Hamiltonian) with the result that angular momentum is 
conserved. 
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ENERGETIC STABILITY 

To be concrete, we will consider the problem of a 
doubly-quantized vortex state. We start by examining 
the static stability of this state within the subspace of 
states of the lowest Landau level. We will generalize 
our arguments below. The corresponding nodeless eigen- 
functions of the (two-dimensional) single-particle Hamil- 
tonian, hi, will be denoted as with m the angular 
momentum; their energy eigenvalues are e^. 

We now consider the state 

4> = CQ(j)Q + C202 -I- C^4>A, (2) 

which describes a doubly-quantized vortex state when 
C2 — 1 and Co = C4 = 0. Since we wish to consider 
the effect of small admixtures of states with m = and 
m = 4, we expand the energy to second order in cq and C4 
to obtain the expectation value of the energy per particle 

E ^ E2 + CqCoMoo + cXciMii -f Afo4(cQC4 -f C0C4), (3) 

where 

E2 = £2 + 2-9^2222 

^^00 = eo ^ £2 + 5(2-^0202 ^ -^2222) 
M44 ~ £4 — 62-1- 3(2/2424 — ^2222) 
Mo4 = 5/0422 = g/2204, (4) 

with Inmik = / (l)n(l>mMkdr and g = NUq. The ele- 
ments, Mmm', are real due to the Hermiticity of H. 
The energy E is then equal to 

/^ = /^2 + I Co P A/00 I C4 P M44 -f 2 1 Co 1 1 C4 1 Mo4 cos (6I0 6*4 ) , 

where Cm = Icmle**"*. The absence of terms linear in cq 
and C4 is a consequence of the spherical symmetry of H 
and implies that C2 = 1 is a solution to the variational 
problem in this restricted space. For fixed |co| and |c4|, it 
is elementary that E has extrema for cos(6'o -I- 64) = ±1- 
In general, the linearly independent curvatures of this 
quadratic energy surface are obtained as the eigenvalues 
of the real symmetric matrix 

, . f Moo Moi \ 
= \ M40 M44 ) 

with M40 = mI^. 

In this problem, it is evidently of physical interest to 
consider the stability of the doubly quantized vortex sub- 
ject to the constraint that angular momentum is con- 
served. This constraint will be satisfied if I Co I = |c4| = |c|, 
and the energy E then becomes 

E = E2 + 2(7|cp(/o202 + /2424 — /2222 ± /o422)- (7) 



The doubly-quantized vortex will be an energy minimum 
if the coefhcient of |cp is positive. In other words, the 
doubly-quantized vortex will be energetically stable if 

/0202 + /2424 — /2222 > /o422- (8) 

This condition is not satisfied when only lowest Landau 
levels are retained 0, and a doubly-quantized vortex 
state is energetically unstable within this approximation. 

Several comments are in order before turning to the dy- 
namical behavior of the same (truncated) problem. We 
have tested the stability of the state with respect to the 
admixture of states (j)o and 4>4- A complete test of sta- 
bility requires the investigation of the admixture of (pm 
for all TO. Fortunately, the spherical symmetry of H en- 
sures that the generalized quadratic form of Eq. ^ will 
only mix states with to = 2 ± fc. With the inclusion of 
additional values of k, the matrix Mg becomes block di- 
agonal and leads to a sequence of equations for each k 
completely analogous to those considered above. (The 
k = 2 choice considered here is known to be the most un- 
stable case jl].) The proble m of finding the eigenvalues 
of the Hermitean matrix, M^, is subject to familiar vari- 
ational arguments. Expansion of the dimension of this 
matrix, either by including more values of k or by relax- 
ing the restriction to the lowest Landau level, can only 
decrease the smallest eigenvalue of Mg. Thus, the fact 
that a given state is energetically unstable for a given 
choice of the finite space used to construct Mg is conclu- 
sive proof of instability. Finally, we note that changing 
the state under investigation, e.g., through the inclusion 
of higher Landau levels, leads to fundamental changes in 
Mg. Under such circumstances, simple variational argu- 
ments cannot tell us whether the improved state will be 
more or less stable. 



DYNAMIC STABILITY 

As mentioned above, dynamic stability probes the tem- 
poral behavior of the system. We again start with the 
wave function of Eq. Q and construct the Lagrangian to 
second order in the small parameters Cq and C4. Varia- 
tions with respect to the coefficients Cq and C4 then lead 
to the equations 

ico = cqCo + 2.gcoc5c2/o202 + 5C4C2C2/0422 

-icl = e^cl + 25C4C2C2/2424 + 5C0C2C2^0422 

ZC2 = e2C2 -I-5C2C2C2/2222, (9) 

where the dot denotes a time derivative. We can solve 
these equations with the ansatz co{t) — co(0)e~*'^'^"'"'^^* 
and cl{t) = c|(0)e^(^-")*. (The time dependence of 
the solution, (j)2, is given by the factor e~*'^* with fi = 
£2 + 5-^2222-) The resulting Bogoliubov equations can be 
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written in the form 

f eo - M + 25/0202 3^0422 ^ '^0 ^ ^ 

V -5^422 -£4 + M - 2g/2424 / V '^4 / 

..(:^).(io) 

This equation can also be written as 

where Ms is the matrix appearing in the static stabil- 
ity problem. The dynamic problem is thus governed 
by the eigenvalues of the non-Hermitean Eq. Hl()(l . and 
these eigenvalues can be complex. Since Mg is real Her- 
mitean, the dynamic eigenvalues are either real or come 
in conjugate pairs. For each eigenfunction, Vdj with a 
complex eigenvalue, Wd, i^d will also be an eigenfunction 
with eigenvalue w^. As a result, the roots move along 
the real axis, touch (i.e., become degenerate) and then 
move off in the complex plane. Evidently, the existence 
of complex eigenvalues indicates exponential divergence 
and dynamic instability. 

From Eqs. pO|l or we see that the Bogoliubov 

eigenvalues are real and the state under investigation is 
dynamically stable provided that /0202 + ^2424 — ^2222 > 
/o422- This condition is identical to the condition for 
static stability found in Eq. JSjl. As we show below in 
greater generality, static stability always implies dynamic 
stability, but dynamic stability can occur even in the 
presence of static instability. 

GENERAL RESULTS 

At this point it is useful to generalize our static and 
dynamic formalisms to allow for the inclusion of an arbi- 
trary number of Landau levels in the description of both 
the doubly quantized vortex and the potentially unstable 
states with m = and m = 4. The matrix appropriate 
for the static problem now has the Hermitean block form 

M - ( ^^04 ^ (12) 
'""^ - M40 M44 ) ' ^'^> 

with Ml = MJ = Ms- Aside from its dimension, the 
only change in the construction of Ms lies in the replace- 
ment of the lowest Landau state, (/)2, by a superposition 
of Landau states, ■02, which satisfies the obvious Euler 
equation 

^oV'2 + ff|V'2p'02 = ^"02- (13) 

Here, /i is a Lagrangian multiplier introduced to ensure 
that ■02 is normalized to unity. It is understood that the 
replacement 4>2 ^2 is to be made as appropriate in the 



integrals, Inmik, contributing to Mg. As before, positive 
eigenvalues of Ms imply energetic stability. 

The dynamic Bogoliubov equations again assume the 
form of Eq. and can be written as 

where Ms is again the real Hermitean matrix governing 
energetic stability. The quantities cq and C4 now repre- 
sent column vectors; their dimensions need not be equal. 
The reality of all eigenvalues, cudi is again an indication 
of dynamical stability. 

As we have seen, H is Hermitean and spherically sym- 
metric. The time evolution of an arbitrary wave func- 
tion, ijj, is governed by the full time-dependent Gross- 
Pitaevskii equation 

hoi: + g\^P\^i> = ih^. (15) 

It thus comes as no surprise that energy and angular mo- 
mentum are, in general, constants of the motion. On the 
other hand, given that the (linearized) Bogoliubov equa- 
tions do not even conserve probability, it is not obvious 
that energy and momentum are conserved when these 
equations are used to describe the temporal evolution 
of the system. This point merits some attention. Con- 
sider a dynamic eigenvector of the Bogoliubov equations, 
ipd, and its (possibly complex) eigenvalue, ujd- Given the 
Hermiticity of Ms, standard arguments reveal that 

iLOd-co:){M'yHd) =0. (16) 

This equation is evidently trivial when uJd is real, but it 
shows that {'4'd\<^\''Pd) = when uJd is complex. This tells 
us that the probabilities of finding m = is equal to that 
for TO = 4 for all times. In this sense, angular momen- 
tum is rigorously conserved in spite of the approximations 
leading to the Bogoliubov equations. 

The extension of this familiar conservation law pre- 
sumes that the trial state ■02 is deformed by the inclusion 
(with arbitrary amplitude) of a single Bogoliubov eigen- 
vector with complex eigenvalue. Angular momentum is 
not in general conserved with arbitrary deformations of 
■02 involving either single Bogoliubov eigenvectors with 
real eigenvalues or superpositions of Bogoliubov eigen- 
vectors. Since the primary rationale for studying the 
Bogoliubov equations at all is to determine the existence 
or non-existence of complex eigenvalues, this point is of 
interest in spite of such caveats. Further, it suggests that 
we are most likely to reveal relations between the prob- 
lems of static and dynamic stability if we consider the 
question of static stability subject to the constraint of 
constant angular momentum. Indeed, this constraint was 
imposed trivially in Eq. ((JJ) above, where it was crucial in 
establishing the identity of static and dynamic stability 
criteria in the simple case of single Landau levels. 
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For the more general problem considered here, it is 
easiest to proceed by introducing real Lagrange multipli- 
ers, A and 77, in order to impose the constraints of over- 
all normalization and constant angular momentum, re- 
spectively, in the static problem. The constrained static 
eigenvalue problem then reads 

M,^, - A J ° ^ V^,. (17) 

In practice, the real parameter ry is adjusted so that 
the resulting static eigenfunctions ips satisfy {ipslclips) — 
0. The desired extrema of (M^) then follow as Cos = 
(■03 1 -Ms I '0s)- The similarity between the constrained 
static problem of Eq. I|17|) and the dynamic eigenvalue 
problem of Eq. H14|l is now obvious. We shall exploit this 
connection below. 



CONNECTION BETWEEN STATIC AND 
DYNAMIC STABILITY 

Given the fact that ('0d|o'|'0d) — 0, we see that 
the inner product of (■0^1 with Eq. H14|) also yields 
(V'd|-Ms|0d) = 0. If ^"2 is deformed by the inclusion (with 
arbitrary amplitude) of a single Bogoliubov eigenvector 
with complex eigenvalue, the energy is rigorously con- 
served and identical to that obtained for the pure state 
"02 ■ We can use any 0^ with complex eigenvalue as a 
trial function to provide a variational upper bound of 
zero for the lowest constrained static eigenvalue, w^. Said 
somewhat more simply, it is elementary that the small- 
est ujs cannot be positive if there exists a state 0(j such 
that ('0d|Ms|0;d) = 0. (Henceforth, we will consider only 
the constrained static problem and will refer to it sim- 
ply as the "static problem".) Thus, if the Bogoliubov 
problem has complex eigenvalues, the lowest eigenvalue 
of the static problem must be negative (or zero). Con- 
versely, if all lbs > 0, none of the Bogoliubov eigenval- 
ues can be complex. In other words, dynamic instability 
guarantees static instability, and static stability guaran- 
tees dynamic stability. This establishes the first rigorous 
connection between the two stability problems. While 
this result may appear obvious, a similar argument shows 
that it is equally impossible to find a state 0^ such that 
(0d|Ms|0d) = if all < 0. This yields the more sur- 
prising result that complete static instability also leads to 
dynamic stability. 

In order to find additional ties between the static and 
dynamic stability problems, it is useful to follow the tra- 
jectory of two Bogoliubov eigenvalues (as a function of 
the coupling constant) as they evolve from distinct real 
values to a conjugate pair. To see this, it is sufficient to 
consider the two-dimensional space of ("00, 0) and (0, tpi). 
(Here, 0o and 04 are of arbitrary dimension and normal- 



ized.) We thus write 

"■-an)- w 

with all quantities real. The Bogoliubov problem of 
Eq. H14I) is readily solved with the following results. For 
\g\ < gc — \{Eo + E4)/2V\, the uJd are real and distinct. 
The eigenvectors are real and do not satisfy the constraint 
of constant angular momentum. For \g\ — g^ the eigen- 
values are degenerate with utd — {Eq — E4) /2. The eigen- 
vectors are identical and given as (■0o,-'04)/\/2. (Recall 
that the problem is non-Hermitean.) For g > gc, the 
eigenvalues form a conjugate pair. The eigenvectors also 
form a conjugate pair with a non-trivial phase, and the 
magnitudes of their V'o and '04 components are equal (i.e., 
they conserve angular momentum). 

We can also solve the static problem using Eq. H17|l in 
the same truncated basis. The values of rj required to im- 
pose the constraint of conserved angular momentum are 
found to be 77 = {E4 ± gV)/{Eo + gV). The correspond- 
ing static eigenvalues are uis = (Eq + i?4)/2 ± gV with 
V's = (V'o,±V'4)/v^j as expected. For \g\ < gc, the static 
eigenvalues are positive. For \g\ > gc, one of the static 
eigenvalues is negative. For \g\ = gc, one static eigen- 
value is precisely zero with = ipd = (V'o, — '04)/\/2- 
The condition for a static eigenvalue of zero is the same 
as the condition for degenerate dynamic eigenvalues. 

We also see that the static and dynamic eigenvectors 
are identical at this point. This should come as no sur- 
prise. At \g\ — gc, the Lagrange multiplier, rj, is — 1. 
Equation (|17|l . which determines ■0s, becomes identical 
to the BogoHubov equation, Eq. ifT^ . and identical solu- 
tions must result. This result is general. While the form 
of these two equations is similar, the fact that Mg is a 
real symmetric matrix forces 4's and uJs to be real (up to 
a trivial phase). When 7^ w^, it is clear that the asso- 
ciated dynamic eigenvectors, ipd and -0^ cannot be chosen 
equal. Thus, the phase of tjjd is non-trivial, and tpd can- 
not be the desired solution to Eq. l(T7|) . (The case of uJd 
real can be dismissed summarily, since the amplitudes of 
TO = and 4 states are not equal.) When two Bogoliubov 
roots are degenerate, however, ipd satisfies both the con- 
straints of angular momentum conservation and reality. 
Under such conditions, -08 = ipd, and Ug = 0. In other 
words, the number of conjugate pairs of complex Bogoli- 
ubov eigenvalues changes by one when a static eigenvalue 
passes through zero. 

It might be thought that the number of conjugate 
pairs of dynamic eigenvalues was equal to the number 
of negative (Dg and that a stronger statement was thus 
possible. As we will show below, this is not the case. 
The present statement can nevertheless give us a use- 
ful corollary. Start from a manifestly stable choice of 
■02 for which all uJs > and dynamical stability is en- 
sured. Vary some external parameter (e.g., the coupling 
constant, g). As demonstrated above, a single conjugate 
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pair of Bogoliubov eigenvalues will appear when the first 
(Ds passes through zero, and the system becomes dynam- 
ically unstable. With every subsequent passage of an 
Qg through zero, the number of conjugate pairs changes 
(i.e., increases or decreases) by one. It is then clear that 
an odd number of iDg must correspond to an odd num- 
ber of conjugate Bogoliubov pairs and, hence, dynamic 
instability. 

DYNAMIC STABILITY IN THE ABSENCE OF 
STATIC STABILITY 

As mentioned earlier, it is possible for the system to be 
dynamically stable even in the presence of static instabil- 
ity. This effect is more subtle but is also useful in tight- 
ening the connection between the problems of static and 
dynamic stability. A simple calculation will be helpful. 
Start with two solutions to the constrained static prob- 
lem of Eg. (|17|l . |i) and |2). The corresponding static 
eigenvalues are (Dsi, which is assumed to be negative, 
and U)s2, which can pass through zero. Approximate the 
Bogoliubov equations by truncating the basis to include 
these two states and the related states |3) — a\l) and 
1 4) = a\2). (This calculation is exact when the dimen- 
sion of Ms is equal to four.) Since these states are not 
necessarily orthogonal, the Bogoliubov equations assume 
the form 

Mg^d = uJd^^d, (19) 

with {Ms)ij = (ilMslj) and (i;)y = Both matri- 

ces are real symmetric. The first two diagonal elements 
of Ms are the constrained static eigenvalues, u>si and lJs2- 
The remaining two diagonal elements, uJsS and Lds4:, will 
be assumed to be positive. The matrix, E, is intimately 
related to the overlap matrix with elements Since 
the basis states conserve angular momentum, the diag- 
onal elements of E vanish. Given the origin of the |T) 
and 1 2) as constrained extrema of (Ms), it is clear that 
the elements of Ms and E are not independent. We shall 
ignore this fact in the following qualitative arguments. 

First, consider the case when these four states are max- 
imally linearly dependent. State |1) has been assumed to 
be distinct from |2) and is explicitly orthogonal to |3). 
Let us assume it is identical to |4). This immediately im- 
plies that 1 3) is identical to |2). The approximate Bogoli- 
ubov equation is thus reduced to a 2 x 2 matrix equation 
in the space of |1) and |2) and assumes the form 

WW)\ ^si {l\Ms\-2))'^' = '''^'- 

(20) 

The eigenvalues of this problem are 

= 7T^ ((1|^^«|2) ± V^^:^) , (21) 



and the corresponding (unnormalized) eigenvectors are 

V'd- v^|i) + v5;t|2). (22) 

When LJs2 > 0, the Bogoliubov eigenvalues form a con- 
jugate pair. Further, Eq. H22|l shows that V'd acquires a 
non-trivial phase and that (ipdWli^d) R-e[(a}*ia}s2)"^^^] 
is zero. The system is dynamically unstable. When 
LUs2 — 0, the Bogoliubov eigenvalues are degenerate and 
real. In this case, ipd is equal to |2), which is real and 
conserves angular momentum. These results are all con- 
sistent with those found above. When ujs2 < 0, however, 
the Bogoliubov eigenvalues are real and distinct. The 
corresponding ipd can now be chosen real, and angular 
momentum is no longer conserved. The presence of two 
negative static minima of (Ms) is thus capable of creat- 
ing dynamic stability in spite of manifest static instabil- 
ity. This argument is actually more detailed than neces- 
sary. We have already seen that the number of conjugate 
pairs of Bogoliubov eigenvalues necessarily changes by 
one every time a static eigenvalue passes through zero. 
Linear dependence has reduced the present problem to 
one spanned by a two-dimensional space. The fact that 
ujsi < ensures that there are initially two complex 
eigenvalues in the space. We have seen in general that 
the number of complex eigenvalues must change by two 
when (Ds2 crosses zero. Since all eigenvalues of this 2x2 
problem are already complex, the only possibility is thus 
that the number of complex eigenvalues is reduced to 
zero with dynamic stability as a consequence, as we have 
seen. States |1) and |2) are both essential to this process. 

Now consider the case where the \i) are all mutually 
orthogonal and maximally linearly independent. Parti- 
tion the matrix into 2x2 blocks spanned by the states 
|1) and 1 3) and the states |2) and |4), respectively. Each 
of the diagonal sub-blocks of Md now has a form familiar 
from Eq. (|20|l . The diagonal elements are equal in each 
case. The products of off-diagonal elements are WsitDsS 
and (Z's2Ws4, respectively. If the off-diagonal blocks of 
this matrix are sufficiently small, the Bogoliubov eigen- 
values will be given by the eigenvalues of these two 2x2 
matrices. Since Wsiu)s3 < by assumption, this block 
produces the original conjugate pair of dynamical eigen- 
values independent of the properties of |2). As ujs2 passes 
through zero, u)s2<^s4 becomes negative, and an additional 
conjugate pair of eigenvalues appears independent of the 
properties of |1). This behaviour will persist whenever 
the off-diagonal blocks of Md are sufficiently small. 

We thus see that there are two possible outcomes 
when a second static extremum becomes negative. If 
the corresponding states are "strongly" coupled with 
|(1|ct|2)| — |(1|4)| Ri 1, the initial dynamic instability 
will be eliminated. If these states are "weakly" coupled 
with |(T|tT|2)| « 0, a second unstable dynamic mode will 
appear. There is another and potentially fruitful way to 
distinguish these alternatives. Consider the evolution of 
the closed surfaces with (Ms) = which bound domains 
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within which {Ms) < 0. [It is assumed that (Ms) is ex- 
plored with real trial vectors subject to the constraint 
of constant angular momentum. Hence, these vectors lie 
on a (hyper)torus.] Every negative extremum of (Ms) is 
evidently surrounded by such a surface, and every such 
surface must contain at least one negative extremum. If 
there is one negative extremum, there is one (Ms) = 
surface. There are two possible ways for the topology of 
these surfaces to change when a second static extremum 
becomes negative. A second (Mg) = surface, not con- 
nected to the first, can appear. In this case, the new neg- 
ative extremum must be a local minimum. Alternatively, 
the original surface can spread, exploiting the fact that 
the manifold is periodic, and touch itself. The point of 
contact necessarily represents a new extremum of {Mg), 
and a few moments thought reveals that this extremum 
is necessarily a saddle point. (When there are more than 
two negative static modes, independent (Mg) = sur- 
faces can merge in a similar manner.) Randomly drawn 
numerical examples suggest that the appearance of a new 
(Ms) = surface corresponds to the case of weak cou- 
pling discussed above, and a new unstable dynamic mode 
emerges. The merger of an (Mg) = surface with itself 
corresponds to the case of strong coupling, and the orig- 
inal unstable dynamic mode disappears. Although the 
evidence is admittedly slim, it is tempting to guess more 
generally that each closed surface with {Mg) = pro- 
duces zero unstable dynamic modes when it contains an 
even number of negative extrema of (Mg) and one un- 
stable dynamic mode when it contains an odd number of 
such extrema. 

As we have shown, the reality and Hermiticity of Mg 
and the non- Hermiticity introduced by the specific form 
of a are not sufficient to provide a unique answer to 
the question of dynamic stability in the presence of an 
even number of unstable static modes. (As noted above, 
an odd number of negative static modes always implies 
dynamic instability.) The question of what constitutes 
"strong" or "weak" coupling would thus seem to require 
more details regarding the physical system in question. 
For this reason, we now turn to the question of the sta- 
bility of vortex solutions to the physically relevant Gross- 
Pitaevskii equation. 

NUMERICAL RESULTS 

We now study the static and dynamic stability of a 
doubly quantized vortex numerically. The single-particle 
Hamiltonian, ho = -h'^V'^ / {2M) + Mcu'^r'^ /2, represents 
a harmonic oscillator potential of strength uj. In an an- 
harmonic potential, a doubly quantized vortex is energet- 
ically stable for weak couplings, but in a purely harmonic 
potential this is known not to be the case 0, . 

In our simulations Eq. H13|) was first solved to find the 
doubly quantized vortex. This is equivalent to finding a 
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FIG. 1: Eigenvalues of the static and dynamic stability ma- 
trices as functions of the coupling constant g, computed by 
means of exact diagonalization in a truncated basis. Solid 
lines represent the imaginary part of the eigenvalues of the 
dynamic stability matrix, lm{ujd)- Dots represent the lowest 
eigenvalues of the static stability matrix, Cus- 

stationary solution to the Gross-Pitaevskii equation sub- 
ject to the constraint that the system has a An phase 
singularity at the origin. This minimization was carried 
out within a truncated basis composed of the Nc lowest 
radially excited states and is valid for small values of g 
such that g/{2hijj) < Nc- The static and dynamic eigen- 
values were then calculated according to Eqs. H14|l and 
(|17|l in a basis consisting of N^. radially excited states for 
both the m — Q and m = 4 components. The computa- 
tion of the eigenvalues was carried out in Matlab. The 
results of this calculation are shown in Fig.^ Results 
are shown for Nc = ~ 8, and we have checked con- 
vergence with respect to both Nc and N^ for the range of 
couplings displayed in the figure. The dynamic eigenval- 
ues LOd coincide with those previously found numerically 
inRefs. null. 

The figure clearly shows the correlation between the 
"window" structure of the complex frequencies and the 
eigenvalues, ojg, of the static problem in agreement with 
our general arguments above. For small g such that 
g <^ Tiuj, the properties of the system follow from the 
perturbative analysis of the lowest Landau level as de- 
scribed above: The system possesses one negative static 
eigenvalue and therefore one pair of complex dynamic 
eigenvalues. As the coupling increases, a second static 
eigenvalue crosses zero at g/huj « 6, and the system be- 
comes dynamically stable despite its manifest static in- 
stability. As a third static eigenvalue crosses zero, one 
pair of dynamic eigenvalues again become complex but 
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become real again when a fourth static eigenvalue be- 
comes negative. The numerical findings indicate that as 
g is further increased, a succession of such alternations 
occurs with the result that the system is found to be dy- 
namically stable (unstable) when there are an even (odd) 
number of negative static eigenvalues. 

For the range of couplings, g, studied here, there is 
never more than one complex pair of dynamic eigenval- 
ues. It is also seen that the system is always statically 
unstable. In the limit of large g, well-known results for 
vortices in an infinite system apply In that limit 
the energy of two vortices increases monotonically with 
decreasing separation, thereby implying that at least 
one eigenvalue of the static stability matrix is negative. 
Hence it is reasonable to conclude that the system is stat- 
ically unstable for all values of g. We thus conclude that 
the "window" structure of alternating dynamically stable 
and unstable regions arises from the non-trivial interplay 
between negative modes of the static stability matrix. 

CONCLUSIONS 

The purpose of this paper has been to point out the 
intimate connections between the static and dynamic sta- 
bility of solutions to the Gross-Pitaevskii equations. We 
have shown that (suitably constrained) static stability 
necessarily implies dynamic stability, that the number of 
complex conjugate pairs of dynamic eigenvalues changes 
by one every time a constrained static eigenvalue passes 
through zero, and that an odd number of negative static 
eigenvalues thus implies dynamic instability. Numerical 
investigations revealed that the doubly-quantized vortex 
solution to the Gross-Pitaevskii equation is statically un- 



stable over the full range of coupling constants explored 
but displays windows of dynamic stability. The general 
nature of the arguments presented here suggests that sim- 
ilar connections between the problems of static and dy- 
namic stability are likely to be wide-spread. Further, we 
believe that additional insight obtained by studying both 
stability problems is likely to be well worth the minimal 
additional effort required. 
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